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We study the impact of an additional U{1)' gauge symmetry with flavour-dependent charges 
for quarks and leptons on the LHC flavour anomalies observed in —>■ K*pL^p,~, R{K) = —> 

KpNpL~ /B — > Ke'^e~ , and h ^ fir. In its minimal version with two scalar doublets, the resulting 
model naturally explains the deviations from the Standard Model observed in B —>■ K* and 
R{K). The CMS access in /i —>■ /rr can be explained by introducing a third scalar doublet, which 
gives rise to a prediction for r —>■ 3/i. We investigate constraints from flavour observables and direct 
LHC searches for pp ^ Z' ^ pL~. Our model successfully generates the measured fermion¬ 
mixing matrices and does not require vector-like fermions, unlike previous attempts to explain these 
anomalies. 


I. INTRODUCTION 

The discovery of a scalar particle at the LHC mm with 
properties close to its theoretical prediction within the 
Standard Model (SM) marks its completion as a descrip¬ 
tion of particle physics. While direct searches for physics 
beyond the SM were negative at the first LHC run, there 
are some interesting indirect hints for new physics effects 
in the flavour sector, namely in the decays of B mesons - 
B K*pl+pl- and R{K) = B ^ Kp+pL-/B ATe+e" 
- and in the decay of the Brout-Englert-Higgs boson 
h —)■ pLT. 

Specifically, the deviations from the SM found by 
LHCb [ 5 ] in the decay B K* pT arise mainly 
in an angular observable called [J], with a signif¬ 
icance of 2-3 a depending on the assumptions for the 
hadronic uncertainties m\- This effect can be ex¬ 
plained in a model-independent effective-field-theory ap¬ 
proach by a fairly large contribution of the operator 
Cg^is'faPLb){ji'y°‘p) [SHIO].^ LHCb further observed 
lepton-non-universality in the H-meson decays [T^] 

± 0.036, (1) 

which deviates from the SM prediction = 1.0003 ± 
0.0001 [13] by 2.6 CT. A possible explanation comes again 
in the form of a non-zero new-physics contribution to 
Cg^ - of the same magnitude as the one required by 
B -X K*p+p- li nn - as long as the analogous con¬ 
tribution to the corresponding operator with electrons, 
(71% is small [ISHIH]. 

CMS recently presented the results of a search for 
a lepton-flavour violating (LFV) decay mode of the 
125 CeV scalar h —>■ pr, with a preferred value m 

Br[h ^ pr] = (0.84%°;39) % , (2) 

updating an earlier preliminary result [20] . Since this 
decay is forbidden in the SM it corresponds to a 2.4 cr 


^ According to Ref. m, also underestimated charm effects could 
explain the deviations from the SM. 


deviation. This is particularly exciting because it hints 
at LFV in the charged-lepton sector, whereas we have 
so far only observed LFV in the neutrino sector through 
oscillations. Seeing as the simplest SM extensions that 
can account for neutrino masses and mixing would not 
lead to observable h ^ pr rates, the confirmation of this 
decay would have a huge impact on our understanding of 
lepton flavour. In particular, it would imply potentially 
measurable rates for LFV processes such as t ^ 3p or 
T ^ p'f Models aiming to accommodate or ex¬ 

plain Fq. ([^ rely on an extended scalar sector [30U35] or 
non-renormalizable effective operators [31IMll37j. 

An explanation for h ^ pr typically requires addi¬ 
tional scalars, while an explanation for B —> K*p'^p~ 
requires additional Z' vector bosons (or leptoquarks 
[ 111131131 ]) to generate the current-current interac¬ 
tion {sjaPLb){pj°‘p) [4nH43] . If the Z' couples non¬ 
uni versally to leptons, it can also account for R{K) = 
B —>• Kp'^p~/B —^ ATe+e” [33]. In Ref. [ 33 ] we pre¬ 
sented a model that can resolve all three anomalies, by 
combining the model of Ref. [33] (with gauged 17(1) l„-l^ 
and effective Z' couplings to quarks generated by heavy 
vector-quarks) with the one of Ref. [32] (with gauged 
17(1 )l^__l^ broken by a second scalar doublet). The com¬ 
bined resolution of the three flavour anomalies gave in 
particular rise to a prediction for the rate of r —3/i. 

In this article we want to study the possibility that 
the same effect employed in Refs. [321 133 ] in the lep¬ 
ton sector might also be responsible for flavour viola¬ 
tion in the quark sector. Therefore, also the quarks 
must be charged under the new 17(1) gauge group, 
leading to a model with flavour dependent B and L 
charges [351135] . In this way, the introduction of vector¬ 
like quarks - which are somewhat “artificially” charged 
under — Lt - can be avoided. Furthermore, the model 
can explain the smallness of the Vub and I4b elements 
of the Cabibbo-Kobayashi-Maskawa (CKM) mixing ma¬ 
trix. Since the flavour diagonal couplings to quarks are 
not arbitrary anymore, interesting correlations with LHC 
searches arise. 


The outline of the article is as follows: In the next 
section, we will consider the possible charge assignments 


and the symmetry breaking in the quark sector. Sec. HI 
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is devoted to a phenomenological analysis of the effects in 


quark flavour physics and direct LHC searches. Sec. IV 


extends the symmetry breaking to the lepton sector, al¬ 
lowing for a simultaneous explanation of /i —> ^t. Finally 
we conclude in Sec. [V| In App. 1^ we briefly discuss re¬ 
lated horizontal gauge symmetries. 


II. MINIMAL MODEL WITH TWO SCALAR 
DOUBLETS 

Here we study the minimal model with flavour- 
dependent 17(1)' charges which can give rise to the de¬ 
sired effects in H —)• K *and R{K). To generate 
the masses and CKM angles, at least two scalar doublets 
are necessary. 


A. Charge assignment 

Concerning leptons, we are drawn to the abelian sym¬ 
metry 17(1) L^-L^ ■ It is an anomaly-free global symmetry 
within the SM [SDH52] and also a good zeroth-order ap¬ 
proximation for neutrino mixing with a quasi-degenerate 
mass spectrum, predicting a maximal atmospheric but 
a vanishing reactor neutrino mixing angle |53H55) . Fur¬ 
thermore, since the Z' boson does not couple to electrons, 
i.e. Cf® = 0, one naturally obtains an effect of the ap¬ 
propriate size in R{K) once acquires its preferred 
value from B —>■ K*^^. Therefore, we choose the follow¬ 
ing assignment for the charges Q' of the new 17(1)' gauge 
group for the lepton generations: 

Q'(L) = (0, 1,-1). (3) 

Breaking is mandatory for a realistic neutrino 

sector, and such a breaking can also induce charged LFV 
processes [SS], such as h ^ fir and r —>■ 3/x ElEaET]. 
However, we postpone the discussion of the symmetry 
breaking in the charged lepton sector to Sec. [TV] 

Concerning the quark sector, the first two generations 
should have the same charges in order to avoid very large 
effects in K-K or D-D mixing, generated otherwise un¬ 
avoidably due to the breaking of the symmetry neces¬ 
sary to generate the Cabibbo angle of the CKM matrix. 
Furthermore, the first two generations mix much more 
strongly among themselves than with the third genera¬ 
tion, so the latter seems to be somewhat special. If we 
require in addition the absence of anomalies, we arrive at 
the following charge assignment for baryons: 

Q'[B) = {-a,-a,2a). (4) 

We will later study the phenomenological implications 
of different values of a. To reiterate, the 17(1)' gauge 
symmetry we consider is generated by^ 

Q' = [L^ - Lr) - a(Hi + 52 - 2B ^), (5) 


so the charges are 


— - for u, d,c,s, 

(6) 

2 a 


y for t, b , 

( 7 ) 

0 for e, Ve , 

(8) 

1 for /r, , and 

(9) 

— 1 for r, Vt . 

(10) 


Note that the relative coupling strength to quarks and 
leptons (parametrized by a) is a free parameter because 
— Lt and Bi + B 2 — 2 B 3 are independently anomaly 
free, a S Q is nevertheless necessary to avoid massless 
Goldstone bosons (see below). Although not required 
by anomaly cancellation, we also introduce three right- 
handed neutrinos vr with charges Q'ii'n) = (0, 1, — 1) to 
employ a seesaw mechanism for neutrino masses. Other 
horizontal gauge symmetries can be considered following 
the above reasoning, but prove less useful in explaining 
the LHCb anomalies (see App. 0. 


B. Scalar sector 

Two SU{2)l scalar doublets 4'i and '^>2 are introduced 
to generate viable fermion masses and quark mixing in¬ 
corporated in the CKM matrix. They carry the f7(l)' 
charges 

g'(4'i) = -a and Q'(^- 2 ) = 0, (11) 

respectively. '1^2, being uncharged under the 17(1)' sym¬ 
metry, will generate the diagonal entries of the fermion 
mass matrices as well as the Cabibbo angle. The vacuum 
expectation value (VEV) (4>i) will generate the mixing 
between the third and first two quark generations neces¬ 
sary for a viable CKM matrix. 

In addition (at least) two SM singlet scalars $1 and 
$2 with charges 

Q'($i) = l, Q'($2) = -a (12) 

have to be introduced to break the 17(1)' gauge symmetry 
above the electroweak scale and generate the coupling 
'1'|4'2 necessary for the mixing of the doublets. The VEV 
($ 1 ) will break the L,- symmetry in the right-handed 
neutrino mass matrix relevant for the seesaw mechanism, 
which leads to a valid neutrino mixing matrix [32l[59]. <I >2 
generates the term (4>2)'k|'I'2 in the scalar potential that 
leads to the 'I'l VEV as well as mixing among 4 '2 and 
4 - 1 . 

For a general a € K, the above particle content 
leads to a Lagrangian with conserved x 

U{ 1 )bi+B 2 - 2 B 3 symmetry. Both the — Lr symme¬ 
tries and the Bi + B 2 — 2 B^ are anomaly free and could 


Gauge symmetries that couple to Bi + B 2 — 2B^ have also been 
discussed in Ref. m with a focus on effects in the up-quark 


sector and coupled in the lepton sector not to — Lti but to 
the lepton charge Lg + — 2Lr (which is not a good symmetry 

in the lepton sector). 





3 


be gauged, giving rise to two additional neutral gauge 
bosons Z' and Z" which can mix with the SM Z bo¬ 
son [60] , However, we only want to promote the linear 
combination of Eq. ® to a gauge symmetry in order to 
end up with a single W that couples to quarks and leptons 
simultaneously. To remove the orthogonal global 17(1)" 
symmetry, generated by Q" = a(L^ — L^) + {Bi + B 2 — 
2 B 3 ), we have to introduce couplings that connect e.g. <i>i 
and $2 non-trivially. This is only possible for a G Q and 
in general requires the introduction of additional medi¬ 
ator fields. Let us sketch the simplest examples for a, 
using the fact that the phenomenologically interesting 
region will be 0 < a < 1 . 

• For a = 1/2, the potential already allows for a term 

that breaks the orthogonal global 17(1)", so 
the particle content from above is sufficient to avoid 
Goldstone modes; a VEV for ^2 will induce a VEV 
for $ 1 . 

• For a = 1/3, it is the coupling that breaks 

the accidental global 17(1)" symmetry explicitly. 

• For a = 1/4, no dimension-4 operators can be writ¬ 

ten down with the given scalars that would break 
the global 17(1)". Therefore, one has to introduce a 
third singlet 4)3 with Q'( 4 ) 3 ) = 1 / 2 , which couples 
via $§$1 and A $2 VEV will induce a VEV 

for $3, which will induce a VEV for $ 1 . 


• For a = 1/6, we need $3 with Q'( 4 ) 3 ) = 1/3 to 
couple 4 > 34 >i and 4 > 2$3 similar to a = 1/4. 


The above considerations show that one can easily con¬ 
struct models for various values of a G Q. However, the 
details of this procedure will hardly make a difference in 
our discussion of the phenomenological effects of the dou¬ 
blet scalars and the Z' in the following. Hence, we treat 
a G Q as a free parameter but use a = 1/2 and a = 1/3 
as benchmark values. 

Similar to Refs. [32l [33] we consider the limit of small 
mixing between the heavy singlet scalars 4)^ and the two 
lighter doublets 






j = l,2. (13) 


Here, 2 correspond to the CP-even and the 

CP-odd components, respectively. For heavy unmixed 
singlet scalars, we end up with a restricted two-Higgs- 
doublet model (2HDM) potential of the form 

+ Ai 2 |'I'ip|'k 2 p + Ai2|4>|4'2p (14) 

+ (to? 2^2^1 +h-C.^ , 


where 771^2 cx is generated by the coupling 

/i$24'l4'2, which induces a small vacuum expectation 
value for 471 [32|. We define v = \Jv\+v\ ~ 246 CeV 
and tan j3 = V 2 /vi , which is medium to large in the region 
of interest. The neutral CP-even components ipi and 


0 R • • 

’ip 2 mix with an angle a in the usual 2HDM notation 
to give the mass eigenstates 

h = cos a — sin a , (15) 

H = sin a tp^’^ + cos a tp^’^ , (16) 

while the CP-odd components as well as the charged ones 
mix with /?, 


A = cos f 3 ip 2 '^ — sin j3tpi’^, (17) 

H~ = cos p-ipi^ — sin/Sipi ■ (18) 

Note that in the general 2HDM —tt/2 < a < 7r/2 [5T| . 
We will always assume that h corresponds to the 125 CeV 
scalar discovered at the LHC mm- Cauge bosons and 
leptons have standard type-I 2HDM couplings to the 
scalars (see for example Ref. [5T]1. 


C. Quark masses and couplings 

The ?7(l)'-neutral scalar doublet '^2 gives flavour diag¬ 
onal mass terms for quarks and leptons, while 411 couples 
only off diagonally to quarks: 

=-Qf{s,l^i+Yf,^ 2 )u. (19) 

~ Q/(C/i'I'i + V/(4'2) di -I-h.c. (20) 

Here, Q is the left-handed quark doublet, u is the right- 
handed up quark and d the right-handed down quark, 
while i and / label the three generations. We also define 
the doublets Tj = JtT 2 'L*. The Yukawa couplings Y“ and 
of 4-2 are forced by our charge assignment to have the 
form 


= 


^21 ^22 
0 0 



( 21 ) 


(with q = u,d) and hence allow for the generation of 
the Cabibbo mixing connecting the first two generations, 
while the third generation is decoupled. The couplings 
are given by 

/o oo\ (Q 0 Ub\ 

r = 0 0 0, e'" = 0 0 (22) 

V etc 0 / \ 0 0 0 / 


and lead to the small mixing of the first two generations 
with the top and bottom quarks after electroweak sym¬ 
metry breaking. The quark-mass matrices in the inter¬ 
action eigenbasis are then given by 

^ (sin/3Y“ -h cos PC) = ULm^Ul. , (23) 

^ ^ y ,03 , (24) 

related to the diagonal mass matrices in the physical basis 
= diag(TO„, rric, mt ), (25) 



4 


= diag(md, ms, mb) (26) The Lagrangian describing the couplings of quarks to 

the physical scalar fields is given by 

by the unitary matrices and The CKM ma¬ 

trix is then given by the misalignment of the left-handed 
up and down quark rotations as 

= uIDl . (27) 


£ D — u 


— u 


COSO 


m^ — 


cos(a — /3) 
•\/2sin/3 
sin(a — /3) 




C ) PRuA + id 

I Pl + 


(28) 


where we omitted the hermitian conjugated terms and 
flavour indices. In addition we defined the non-diagonal 
coupling matrices 

r = uICUr , = dI^^^Dr . (29) 

The terms proportional to m^ correspond to the usual 
tj^pe-I 2HDM-like couplings, while the terms involving 
induce flavour violation and will be regarded as small 
perturbations with respect to the type-I structure. 

Since and correspond to right-handed rotations, 
the CKM matrix is (to a good approximation) given by 
V ~ Dl , while the mixing angles in Dr generated by ^^3 
and ^23 &re suppressed to those in by md^s/^b, and 
so Dr ~ I. a perturbative diagonalization of the quark 
mass matrices gives 


V2 mb V2 mb 

^db — Q ^ub 1 ss6 — o 5 


cos/? V 


cos/? V 


(30) 


in particular ^dhl^sb ^ 1- For medium to large values 
of tan/?, our model explains why Vub and Vcb are much 
smaller than the Cabibbo angle as the contributions of ^sb 
and S^db to the mass matrix are suppressed by tan /?.^ The 
non-diagonal matrix relevant for all the scalar couplings 
is then given by 


I" 




■\/2 mb 
cos/? V 


0 0 -v;,Vtb \ 
0 0 -V*Vtb 
\o 0 i-\vtb\y 


(31) 


The dominant flavour off-diagonal coupling is that to sb, 
while that to db is suppressed by VtdiVts — 0.2. Note 


^ Demanding perturbatively small values for ^db sb yields a mild 
upper bound tan/3 < 0(10^) due to the smallness of Vub,cb- 


that the flavour-violating couplings are entirely given in 
terms of the CKM matrix elements. Therefore, flavour 
violation in the down-quark sector induced by the scalars 
only involves the free parameters tan/?, a and the scalar 
masses, making our model highly predictive. 

While the couplings ^23 ^^3 are fixed by the re¬ 

quirement that the measured CKM matrix is generated, 
the couplings ^33 and Q 2 free parameters generating 
top-quark flavor-changing effects not under investigation 
in this analysis. We therefore neglect their effect in the 
following and set = 0. One should keep in mind that 
these couplings can in principle induce decays such as 
t ^ he and t —)■ hu. However, as we will see, the scalar¬ 
mixing angle a must be very small in order not to vi¬ 
olate bounds from Bg-Bg mixing, rendering the h-q-q 
couplings nearly SM-like. Therefore, also the effect in 
t ^ he and t ^ hu is suppressed. 

The Z' couplings to up quarks are vector-like in the 
limit ^“ = 0 we consider in this paper 

g' uy^diag)—a/3, —a/3, 2a/3)u Z'^ . (32) 

Due to the rotation D^ ~ V, the Z' couplings to down 
quarks are chiral and given by 

g' {da^PLddZ'^Vfj^ + da^PRd.Z'^rf) , (33) 

with coupling matrices 

/\Vtd\^-l VtgV^d VbbV^d \ 

V,dVt*g VtbV,% 

V v.dV^b \Vtb?-\) 

/-1/3 0 0 \ 

T'^^-a 0 -1/3 0 . 

\ 0 0 2/3 J 

Note that this reduces to diag(—a/3, —a/3, 2a/3) also for 
the left-handed couplings in the limit in which the CKM 


(34) 

(35) 
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matrix is the unit matrix. As the right-handed couplings 
are to a good approximation flavour diagonal, this will ul¬ 
timately lead to the desired hierarchy \Cg\ ^ ICg] hinted 
at by global fits [BllHllSl]. Note that an opposite U{iy 
charge for 'I'l would require the CKM matrix to be gener¬ 
ated in the up sector, while right-handed down-quark Z' 
couplings would be free parameters, rendering the left- 
handed down-quark Z' couplings flavour conserving. As 
we will see in the next section, one cannot explain R{K) 
and B —>• in such a setup. Also note that the cou¬ 

pling matrix is hermitian with complex off-diagonal 
entries. This will in principle lead to quark-dipole mo¬ 
ments |63j . which are, however, proportional to F^-F^* 
(with i j) and hence vanishingly small in our scenario 
with diagonal F-^. 

It is important to reiterate that our model has, by con¬ 
struction, suppressed flavour violation between the first 
two generations, which drastically softens constraints. 
Furthermore, the “detached” third quark generation is 
motivated by the observed smaller mixing angles. 


D. Lepton masses and couplings 


Since we only consider |a| < 1 in order to avoid strin¬ 
gent constraints from AF = 2 processes, it is not possible 
to couple the scalar doublet 'I'l with t/(l)' charge a to 
leptons. The Dirac mass matrices for the leptons are 
hence generated by 'Fg only and are diagonal in flavour 
space due to the charge assignment of our U{1)' symme¬ 
try. There is hence no charged-lepton flavor violation in 
this simple 2HDM. However, extending the model by a 
third doublet with \Q'\ = 2 (Sec. IV) can again lead to 
LFV and in particular give rise to the decay h —>■ /rr, as 
pointed out in Refs. [32l[33]. For now, we consider the 


2HDM with charged-lepton flavour conservation. 

Neutrino masses arise via the seesaw mechanism [55] 
as rrii, ~ —moTny^rny, with diagonal Dirac mass ma¬ 
trix mu oc (^' 2 ) and the right-handed Majorana neutrino 
mass matrix given by |59j 


Ml ai2($i) ai3(‘hi) 

’niR = I ai2($i) 0 M 2 

.ai3($i) M 2 0 


(36) 


As a consequence of our minimal singlet sector, two tex¬ 
ture zeros in mu remain and propagate to m^, as two 
vanishing minors |59j : (m“ ^)22 = 0 = (m^i) 33 . These 
zeros lead to relations among the neutrino mixing pa¬ 
rameters that allow us to predict the unknown phases 
and masses given the measured neutrino mixing angles 
and the mass-squared differences [551 EOj. In our case we 
predict normal ordering, typically with quasi-degenerate 
masses. Delicate cancellations can, however, occur which 
lead to a normal hierarchy (see Fig. With the best-fit 
values for Amf^ and 9ij from the global fit of Ref. |51| , we 
predict the lightest neutrino mass to be m-i = 0.076 eV, 
and hence = 0.24 eV, as well as | sin^cpl = 0.96, 

and mpp = \{m,^)ee\ = 0.065eV. Note that this predic¬ 
tion for 5cp agrees well with the preferred region from 
the global fit, 5 cp/° = SOb^^yg [64]. Using the 3a ranges 


for mixing angles and Amf^, we obtain lower bounds on 
the neutrino mass observables 

mj > 0.095 eV, meV. (37) 

3 


The sum of neutrino masses is large enough to 

have a cosmological impact m-, in fact, Planck already 
constrains < 0.23eV at 95% C.L. which, 

however, depends strongly on the combination of data 
sets. The effective neutrino mass mpp relevant for neu¬ 
trinoless double-beta decay is also in the reach of future 
experiments (see Ref. m for a recent review). There 
is a strong correlation between the CP phase and the 
atmospheric mixing angle, approximately given by 


<5cp 


90° -f 360°(sin2 623 - 1/2), or 
270° - 360°(sin2 6»23-l/2). 


(38) 


For a more extensive discussion of this neutrino structure, 
we refer to Refs. [MUTOj. Despite the lack of charged lep¬ 
ton flavour violation, our simple 2HDM is hence testable 
in the lepton sector as well. 

Let us come to the scalar-lepton interactions. Since 
dti has no couplings to leptons, the couplings of h, H, A, 
and are given simply by those of a type-I 2HDM, i.e. 


£ D — e 


— le 


cos a 
V sin /? 

V tan (3 
V tan /3 


] eh — e 


sm a 
wsin/3 


1 eH 


756 A 

^ Prg + h-c. 


(39) 


with the Pontecorvo-Maki-Nakagawa-Sakata mixing 
matrix U and = diag(me, m,-). Here, we ignored 
the couplings to the right-handed neutrinos, assuming 
them to be very heavy. For the charged-lepton couplings 
to Z', we have the usual vectorial — Lt couplings 


9 ' {ea^PLejZ'^rt^ + ) , (40) 


with 


T^eL _ Y'eR _ 



(41) 


which are flavour-conserving in this minimal 2HDM 
(lepton-flavour violation will arise in the 3HDM of 


Sec. IV). 


E. Gauge boson sector 

The Z' mass in case of just two singlet scalars (see 
Eq. (HI) is given in the limit of interest ($j) ^ (4'^) by 

m%, ~ 25 '^(( 4 >i) 2 -Ha 2 ($ 2 )^). 


(42) 











6 



0.001 


0.001 


0.01 

^lightest 



mpp [eV] 


FIG. 1: Left: mpp vs. lightest neutrino mass (in blue) as predicted by the two vanishing minors {m ~^)22 = 0 = (m “^)33 in 
combination with the 3 (t values for 9ij and Am^j from Ref. [64]. The diamond marks the best-ht point. The dotted black lines 
show the standard allowed Scr range for normal ordering. The red regions are disfavored by GERDA [65] (conservative limit) 
and Planck |66l I67 |. Right: Correlation of sin^ 623 with mpp due to the two vanishing minors. The yellow region denotes the 
Ict range for 823 , and the diamond marks the best-fit point. 


Ignoring kinetic mixing between U{iy and hypercharge 
U{1)y, a Z-Z' mixing angle 6zz' m is nevertheless in¬ 
duced by the VEV of [5^ . 


g'dzz' — — 


agiv'^ cos^ f3 

2m|,/g'2 


~ -2 X 10"'‘a 


10 


tan (3/ xmz'j g' 


TeV 


(43) 


The gauge eigenstates Z and Z' can then be expressed 
in terms of the mass eigenstates Z^ and Z^ as 


Z^Z^-Bzz'Z2, (44) 

Z' ^Z^ + Qzz'Z^, (45) 


easily satisfied for the parameters under investigation 
(see Eq. (43)). 

The Z also inherits the Z' quark couplings and thus 
contributes to flavor violation in the quark sector. Its 
contribution will be enhanced by its smaller mass but 
suppressed by Qzz'- For TeV scale Z', the mass- 
enhancement cannot overcome the Z-Z' mixing suppres¬ 
sion, so we ignore the flavour-changing Z couplings in the 
following. 

Since the Z-Z' mixing angle is small, we will continue 
to denote the new gauge boson by Z' in the following. As 
we assume the three right-handed neutrinos to be heavier 
than the Z', the invisible branching ratio is set by the 
active neutrinos. 


which couple to 

-£ D (51 ji,M + g'^zz’j'f,) Zy + [gj^ - gi0zz'ji,f^) ^2 . 

(46) 

where gi = g/cw = ejswcw (with s\\r = sin0w etc. for 
the weak mixing angle 9w) and the SM neutral cur¬ 
rent [73] 


Br[Z' 


—)■ mv = 


n.nz' 


1 


(50) 


Eyr[Z'^//] 3 + 4a2 

assuming rrit mz' < 2m,y^. In the same limit, we have 

1 


Bi[Z' —>■ fifj] ~ Br[Z' —tt] ~ 


3 -b 4a2 ’ 


(51) 


jT = E [(^ 3 (/) - slrQ{f))PL - sl,Q{f)PR] f . 

(47) 


while the branching ratio into two up quarks (same for d, 
s, c) is a^/3 times the above, that into top quarks (same 
for bottom) is Aa? times the above. 


The U{\)' current is given in Eqs. (32), (33), and ( [40| . 

In the lepton sector, this Z-Z' mixing leads to small 
vector-coupling shifts of the light Z\ (i.e. the SM Z stud¬ 
ied at LEP) to muons and taus, 

(mm, tt) ~ -1/2 -b 2s^ ± g'dzz'lgi , (48) 

and thus ultimately to lepton-flavour non-universality in 
the Z couplings [56j. Simply demanding that the addi¬ 
tional contribution does not exceed three times the Ict 
error {gyirr) = —0.0366 ± 0.0010) gives 


\g'0zz' \ < Ml 0.003 ~ 2 X 10 


-3 


(49) 


III. FLAVOUR OBSERVABLES AND LHC 
CONSTRAINTS 


We will now investigate the relevant bounds on our 
model, i.e. i3-meson decays, AF = 2 processes, neutrino 
trident production, and direct LHC searches. We will 
not go into details about the standard phenomenology of 
2HDMs of type I (see Refs. [TJ] |75| for a recent evalua¬ 
tion) but rather discuss constraints arising from the de¬ 
viations from the type-I structure due to the additional 
flavour violation. 
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A. B ^ K*n+ij,- and B^Kfi+fi-fB^Ke+e- 


Concerning leptonic B decays, both B —>• 
and R{K) are sensitive to the Wilson coefficients 
and incorporated in the effective Hamiltonian, 


BeS = — 


with 


dGp 

TT 


VtbV, 



cfof + cf ^o'" 


+ h.c.. 


(52) 


^ [srPLb] [h^i] , 

Of = ^ [srPRb] ] , 

Ofo = ^ [^"7^7"^] , 

Of = ^ [s7^Pk6] [^7m7'^] , 


with i G {e,/r,r}. In our model we generate the two 
coefficients Cf = —Cg'^ with 


M pTeVy 

® ctem Of Vl/3/ x^Rz'lg’} 

(54) 

while 

C'l'^ = Cf = = Cf = 0. (55) 


Note that Cf is real in our model but could have either 
sign depending on a. As already noted in Refs. 0162 ], 
Cf < 0 and Cf ^ = 0 gives a good fit to data. Using 
the global fit of Ref. 0 we see that at (1 cr) 2 a level 


- 0.60 (-0.95) > Cf > (-1.65) - 2.00, (56) 

with a best-fit value Cf ~ —1.3. Interestingly, the re¬ 
gions for Cg^ required by R{K) and B -G fie 

approximately in the same region. A Z' explanation of 
R{K) and B -G K*jjL^fjT has also been considered in 
Ref. 0 as an effective model, with similar constraints as 
presented here. 

Putting everything together, we find the 2 a preferred 
region for our model from R{K) and B -G A'*/i+p,“ 


2.3 


(10TeV)2 


< ag'^ < 7.7 


(lOTeV) 


2 ’ 


(57) 


with best-fit value mz' j\fag' = 4.5 TeV. In the following 
we discuss additional constraints on a, g' and mz' that 
are of relevance to our resolution of the LHCb anomalies. 


B. AU = 2 processes 

The Z' and scalars have flavour non-diagonal cou¬ 
plings, generating tree-level contributions to AU = 2 pro¬ 
cesses. Considering for definiteness Bg-Bg mixing {B^,- 
Bd and K-K follow by simple replacements of indices). 


the relevant effective Hamiltonian is conventionally writ¬ 
ten as (see for example Refs. [ZillT]) 


= E C'. O, +^C' O' + h.c., (58) 

1=1 j=i 

with 

Ol = [Sa"f^PLba] [sisYPLbid] , 

O 2 [^ckPl^ck] [^fsPLb^] , 

O 3 = [sc^PLbp] [spPLb ^], (59) 

O 4 = [SaPLba] [spPRbp] , 

O 5 = [SaPLbp] [spPRba] ■ 


Here, a and j3 are color indices, and the primed operators 
are obtained by the exchange L gg R. 

In our model we find the contributions (using the ex¬ 
pression for C 2 given in Ref. [77]) 




dL\ 

b j 


O'= 


E — 


(^L)^ (60) 


while all other Wilson coefficients are zero. It is impor¬ 
tant to note that the Z' contribution to the AU = 2 
processes scales like a^g'"^ while the contribution to 

leptonic decays, such as H —?► Kfig, are proportional to 
ag '"^This makes it possible to evade AO = 2 con¬ 
straints by choosing a ^ 1 , while still accommodating 
the LHCb anomalies. As we will see below, a < 1 is 
actually already sufficient for this purpose. 

The Wilson coefficients Cj enter physical observables, 
i.e. mass differences and CP asymmetries, via the calcu¬ 
lation of matrix elements involving decay constants and 
bag factors calculated with lattice QCD (see for example 
Ref. [7H] for a review of recent lattice values). In ad¬ 
dition, the QCD renormalization group effects must be 
taken into account. For this we use the next-to-leading- 
order equations calculated in Refs. [zaiis]. 

On the experimental side, the central values of Ams^ 
and Ams^ are slightly above the SM prediction, and the 
same is true for ek extracted from K-K mixing. This is 
interesting since our model predicts necessarily construc¬ 
tive interference of the Z' contribution with the SM in 
all three observables. For our numerical values, we use 
the 95% C.L. results of the UTfit collaboration [801182] : 

0.76 < Rb, = AmB,/Am|“ < 1.43, (61) 

0.90 < = Ams./Aml^ < 1.23, (62) 

0.77 < Re^ = SK/ej^ < 1.41. (63) 


Similar results are obtained by the CKMfitter collabora¬ 
tion [83] . 

The strongest constraints come from Bg mixing, as can 
be seen in Fig. [^ but can easily be evaded for a < 1 
even for the large \Cg^\ ^ 1 required to explain the B- 
meson anomalies. Similar but weaker bounds hold for 
the other AF = 2 processes. From Bg-Bg mixing we get 
the approximate 95% C.L. bound 


a^g'^ < 2.6 , 

^ - (10TeV)2 


( 64 ) 
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FIG. 2: Rsq and Rsj^ as a function of C'g'^ for different val¬ 
ues of a taking into account the Z' contribution only. The 
horizontal gray regions are excluded by Ams^ and ek, while 
the vertical ones are excluded by B —>■ and R{K). 

We used = 3TeV for the renormalization group for the 
AF = 2 processes. Note that the dependence on mz' is there¬ 
fore only logarithmic. 


FIG. 3: Allowed regions in the tan l3-a plane assuming that 
Cg^ is reproduced by the Z' contribution within 2 a for a = 
1/3 for mu = 300 GeV and niA = 350 GeV (yellow), jtia = 
300 GeV (red) and mriA = 250 GeV (blue). The gray region is 
excluded by 6 —>■ S 7 |84| . Dashed lines indicate cos(a — /3) = 
±0.4 for reference. 


Note that this bound is obtained for m^/ = 10 TeV and 
only scales approximately like in\i due to the additional 
logarithmic dependence on the mass from the renormal¬ 
ization group. 

The flavour-changing neutral scalar couplings also af¬ 
fect Atub, (and also Atob, and Attik)- Here, the H,A 
contributions decouple as Ijrnj^ while the h contribu¬ 
tions vanishes for a — /3 = 7r/2. For large values of tan /3 
and small a we have approximately 


A h,H,A 
/\Tn^ 

Am^ 


0.12 cos^(a — /3) tan^ /3 

2o/(200GeV)2 

± 0.19 tan^ /3 ( -^- 

V 'm-H 


(200 GeV)^ \ 
m\ ) ■ 
(65) 


D. Direct LHC searches 


Since (unlike in Refs. [531 Si]) we have (potentially 
sizable) Z' couplings to the first-generation quarks, our 
model is constrained by LHG searches for pp ^ Z' ^ 
ppL. With vectorial Z' couplings, universal in the first 
four quark generations, our model is closely related to 
U{1)b-l models |S7|, for which dedicated analyses exist. 
Working in the narrow-width approximation, the relevant 
quantity for collider searches is the Drell-Yan production 
cross section times branching ratio into muons, 

a{pp -A Z')BR[Z' -A pp] 1/(3 ± 4a^) 

a(pp Z'g_^)BR[Z'g_^ -A H 2/13 

(67) 


The allowed regions in the tan (3 — a plane are shown in 
the right plot of Fig. assuming that the Z' generates 
the central value of Cg for a = 1/3. One can see that 
the bounds are weaker if uia < mn and stronger for 
rriA > ttih due to the negative A contribution in Eq. (65). 


C. Neutrino trident production 


The most stringent bound on flavour-diagonal Z' cou¬ 
plings to muons only (i.e. no quark or electron cou¬ 
plings) arises from neutrino trident production (NTP) 
VfiN -A v^NpApi~ gHllS]: 


tJNTP 


1 ± (1 ± 4s^ ± 


92 


o-; 


SM 

NTP 


1 + (1 + 4s^) 


( 66 ) 


Taking only the GGFR data [SS], we find roughly 
/g’ ^ 550 GeV at 95% C.L. for a heavy Z'. This cuts 
slightly into the parameter space allowed by Bg mixing 
and Cg but is only relevant for a <C 1 (see Fig. [^. 


valid for rrit mz' < One can therefore simply 

rescale the B — L limits of ATLAS’s = 8 TeV anal¬ 
ysis of di-muon resonance searches (see auxiliary figures 
of Ref. [88]), resulting in Fig. Gontrary to low-energy 
observables - which only depend on the ratio g'"^ jm^, - 
the LHC probes on-shell Z' . This leads to a complicated 
dependence on mz ’, since the production cross section in¬ 
volves parton-density functions. As a result, the produc¬ 
tion cross section, and therefore the sensitivity, decreases 
strongly once mz' approaches the maximal available en¬ 
ergy. 

For a <C 1 (i.e. the leptophilic case), the Drell-Yan 
production becomes negligible, and our model again re¬ 
sembles the standard L^ — Lr models. In this case, the Z' 
production at colliders goes through pp -A' fifiZ' -A 4/i 
(or 3/i plus missing energy), where the Z' is radiated 
off a final state lepton |S3]. For mz' > mz, the LHC 
constraints are currently weaker than NTP [44] but will 
become competitive with higher luminosities [5DH55] . 

Finally, for Z' masses above the on-shell threshold, one 
can obtain limits from searches for contact interactions. 
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mzlQ' [TeV] 


FIG. 4: Limits on qq ^ Z' ^ fiJL from ATLAS [88] (black, 
allowed region down right) and the 2 a limits on Cgto ac¬ 
commodate B —>■ and B —>■ ij.~/B —>■ Ke^e~ 

(red, allowed regions inside the cone). Solid (dashed) lines are 
for a = 1/2 (a = 1/3). For a = 1/2, the green shaded region 
is allowed (similar for a = 1/3 using the dashed bounds). 


in our case 

a g'^ _ _ 

C D - q'y°‘q , with q G {tt, d, s, c} . (68) 
3 

For positive a, the strongest limit from ATLAS is on the 
operator q^'^Phq providing a 95% C.L. limit 

of [ 55 ] 


mz'/g' > 1.4TeVv/a/(l/3), (69) 


FIG. 5: Limits on mz'lg' vs. a from NTP (gray), Bs-Bs 
mixing (red), and (green). The horizontal lines indicate 
some values of interest: a = 1, 1/2, and 1/3. Not shown are 
LHC limits (see Fig. . 


Ke~^e~ within 2a (Fig. This also implies a lower 
limit on the gauge coupling g' > 0.55 (0.65) for a = 1/2 
(1/3), resulting in a C/(l)' Landau pole below 10^® GeV 
(3 X 10^2 GeV). 


We remark that the dominant flavour violation in the 
b-s sector also induces the decay h -G bs, with branching 
ratio of order 10“^ cos^(a — /3) tan^ /?. While generically 
unobservably small due to the i?s-mixing constraints in 
Eq. (65), it can be large if the A contribution to 
takes just the right value. 


which is weaker than the bounds from Cg (Eq. 


IV. EXTENSION TO THREE SCALAR 
DOUBLETS 


E. Discussion 

The relevant low-energy constraints are collected in 
Fig.i If we want to explain B -G K* gA g, and B -G 
Kg^g~/B —>• iFe+e“ within 2 ct (Ict), we need a < 1.13 
(0.71) to avoid stringent B^-Bg mixing constraints (tak¬ 
ing into account the Z' contribution only). Due to the 
stronger dependence on a, the Bg-mixing constraints are, 
however, unproblematic for smaller values of a, and ac¬ 
tually in agreement with the whole 2 cr range for Cg for 
a < 1/3. Values like a = 1/2 or a = 1/3 and mz' jg' — 2- 
4 TeV can therefore easily lead to the required Cg contri¬ 
bution necessary to explain B -G K*g^g~ and R{K) 
(Fig. §. Note that for these statements we assumed 
tua = ruH, i.e. only took the Z' contribution to Bg-Bg 
mixing into account. However, for rriA < rnn the bounds 
get weakened, while they become stronger for ruA > ttih 
due to the (destructive) constructive interference of the 
H (A) contribution with the Z' and the SM one. 

For a < 1/3 and mz'lg' = O(TeV), direct searches 
at the LHC cut into the mz'-g' parameter space that is 
unconstrained by low-energy processes. We then need 
mz' > 2.55 TeV (2.46 TeV) for a = 1/2 (1/3) if we 
want to explain B -G K*g^g~ and B —)■ Kg'^g~/B -G 


Above we considered a 2HDM with a horizontal 11(1)' 
gauge symmetry that leads to flavour-violating couplings 
of h and Z' to quarks and can successfully explain the 
anomalies in B —)■ K*g'^g~ and R{K). In this section 
we will additionally aim at explaining the tantalizing hint 
for h ^ gr from CMS [T5] (see Eq. ([^), which violates 
— Lt by two units. The signal can be accommodated 
in gauged models by breaking the symmetry 

with a scalar doublet IE's carrying \Q'\ — 2 [ 55 ][ 55 ]. Since 
we cannot set |a| = 2 in our 2HDM from above if we want 
to explain the LHCb anomalies (see Eig. we have to 
introduce a third doublet that carries \Q'\ = 2. Thus, in 
total three scalar doublets. 




{Vj + V' 


j = 1,2,3, (70) 


^ Note that the ATLAS constraints can also be evaded for <C 
TeV with much smaller g' (Fig. [^. For a = 1/2 (1/3), this 
would require m^/ < 300GeV (400GeV) and g' < 0.06 (0.1), 
not necessarily compatible with the approximations used above, 
so we omit a discussion for now. 














10 


with V = y/vf +V 2 +v^ ~ 246 GeV, are introduced to 
generate viable fermion masses and mixing, as well as 
the desired lepton-flavor violation in h —>■ /ir. They have 
the U{iy charges 

Q'(vl/i) = -a, Q'(vI/2)=0, Q'(vI/ 3 ) = -2, (71) 

respectively. We need again (at least) two singlet scalars 
to break the 17(1)' spontaneously above the electroweak 
scale. The same arguments from above apply regarding 
the need for more scalars if a < 1/3 is taken. 

For simplicity we will again assume a negligible mixing 
of the doublets and singlets $, as in Sec. |llj which 
makes sense in the limit (4)^),In this 
case, we effectively end up with a restricted 3HDM at 
low energies with a scalar potential 



+ Ai2|4'ip|4'2p + Ai3|4'ip|4'3p -I- A23|4'2nd>3p 
+ Ai2|4'|'I'2p + Ai3|'I'|'k3p -|- A23|'l'2^3p 
-k + h.c.) . 

(72) 


The mass terms of the last line are of particular im¬ 
portance and are generated by the singlet VEVs ($j); 
note the absence of a term 771 ^ 34 '| 4'3 and quartic terms 
like 'l' 24 'i'I' 34 ' 3 . Even though it may look complicated, 
the above potential (with 16 real parameters) is already 
much simpler than the general 3HDM (with 54 parame¬ 
ters) due to our imposed (and essentially softly broken) 
17(1)' symmetry. 

Compared to the 2HDM from Sec. [llj we now need 
two angles /3 to parametrize the ratios ('I' 2 )/('I'i. 3 ), and 
three angles a to describe the mixing of the three CP- 
even fields h, Hi^ 2 - In addition, also the CP-odd fields 
Ai ^2 and the charged fields 11^2 niix among each other. 

Assuming a CP-conserving potential, i.e. real 777^2 231 
the minimization conditions give 


ml 


X 2 VI 

2 


m 


2 ^ 
12 — 


2 mlvi 

V 2 


m 


2 ^ 

23 — 


2mlv3 

V 2 


(73) 


in the limit V 2 ^ U 1 . 3 . This limit is particularly useful in 
this 3HDM because it allows for approximations leading 


Mi 


I ml + ^{ Xl 2 + Ai2)u^ ^5(Ai2 + Ai2)u 

A2U^ 


V 


to rather simple analytic expressions. It is valid because 
we already know from the above analysis that vi V 2 , 
and from Refs. [32l |33] that V 3 V 2 , i.e. we are in the 
“large tan/?” region in both sectors. 

For the pseudo-scalars, the would-be Goldstone bo¬ 
son eaten by the Z is given by G° = J2j Hi® 

pseudo-scalars orthogonal to G° further mix due to a 
complicated 2x2 mass matrix. We write (Ai, G°, A 2 ) = 
L 7 i(' 03 '\' 02 ’\ with the parametrization 


'-1 


Ua = 


-1 


\ r 


cos 7/23 

sin 7/23 

/ \ - sin 7/23 

cos 7/23 / 


cos 7/13 

sin 7/13 \ 

/ cos 7/12 

-sin 7/12 \ 



sin 7/12 

cos 7/12 

- sin 7/13 

cos 7/13 / 

V 

V 


(74) 


For ?;i _3 <C V 2 , the angles are approximately given by 


bl2 - — , ??23 - — , (75) 

V V 

V1V3 2 m§ -I- (A23 + A23 )u^ 

77l Q — -z:-iic- . 

2 (mf - m§) -k (Ai2 + A12 - A23 - X 23 )v‘^ 

(76) 

Here, 7713 is the most complicated angle, but turns out 
to be small due to Ui 7;3 <C While it is obvious that 
the 12 and 23 mixing is small in the limit ui _3 <C U 2 , it is 
less obvious why the 13 mixing is so small. The reason is 
the absence of a term 777 . 334 '| 4'3 in the potential, which 
implies that all 13 mixing originates from the products 
771773 of VEVs. We will neglect the mixing in the 13 sec¬ 
tor in the following, which drastically simplifies matters. 
In the limit 771^3 <§; 772 , we then define the two beta an¬ 
gles /?i 2,23 - 7’‘/2 - 7/12,23, i-e. tan/?y ~ I/t/^-. /?i 2 then 
corresponds to the /? angle from above (Sec. jn|, while 
P 23 corresponds to the (3 angle from Refs. [351ES] and is 
hence expected to satisfy tan /?23 1 (to be quantified 

below). 

The charged fields 7 /^ have the same mixing pattern to 
linear order in 773 , 3 : (TF/", G+, 77^) ~ [/^(t/:/", 7 / 7 ^, 7 /;^), 
with would-be Goldstone boson G+. Einally, the CP- 
even fields (V’l’^, V' 2 ’^, ' 03 ’^) have a symmetric mass ma¬ 
trix 


- 777 ^/ 77^771 (A 13 -|- Al 3 ) 77 l 773 ^ 

(|(A 23 + A 23)77 - 777^/77^ 773 

ml + \{X23 + X23)v'^ j 


(77) 


The mixing is obviously again suppressed by 771 , 3 / 77 ; in particular the 13 mixing is small. We write (Hi, h, H 2 ) 
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with 


Uh ^ cosa23 

y sin a 23 

(78) 

where we already neglected the small 13 mixing. The 
remaining two scalar mixing angles take the form 


- Sm 023 
cos 023 


C0SO12 SmOi2 

— sinoi 2 cos 012 


vi 2 m\ - (Ai 2 + Ai2)n^ 

V 2m\ — 2 \ 2 v‘^ + (Ai 2 + 3 ^ 12 )^^ ' 

^3 2m§ - (A23 + A 23 )^^^ 

V 2m§ — 2\2v'^ + (A23 + \2z)v‘^ 


(79) 

(80) 


Keeping in mind that the second and third scalar dou¬ 
blets hardly mix in the limit of interest, our 3HDM sim- 
plihes significantly, especially taking into account that 
' 1>3 only couples to leptons and d>i only to quarks. Es¬ 
sentially, our model looks like the 2HDM from Ref. [32] in 
the lepton sector, and like a separate 2HDM in the quark 
sector. This means we can describe the lepton sector as 
a type-1-like 2HDM with an angle ,023 and 023 , and the 
quark sector as a type-I-like 2HDM with angle P 12 and 
012 . The scalar doublet ^'2 is mostly SM-like and will 
inherit off-diagonal couplings to /ir from d >3 and to tu, 
tc, db, sb from 'hi. 

In this limit, our 3HDM is parametrized by the masses 
(m/i, WAi.a)) the vacuum-angles /3i2,23, 

and the mixing angles 012,23 for the CP-even scalars. 


A. Quark masses and couplings 

For the quark masses and couplings, nothing changes 
compared to the 2HDM above; we simply replace tan /3 —>■ 
tan/ 3 i 2 — v/vi and o —>■ 012 , keeping in mind that this 
only works in the limits of large tan jdij and small Oy. 
In particular, we can again easily resolve the anomalies 
in R —>• K*p,'^pL~ and B —)• Kp'^p,~/B —)• Ke+e“ using 
the flavour violating Z' couplings as discussed above. As 
we will see below, the additional constraints from the 
lepton sector (compared to the 2HDM discussed before) 
do not interfere with this solution. In fact, the additional 
resolution of /i —>■ /ir makes possible a prediction for 
T —>■ 3/i, depending on Cg. 


B. Lepton masses and couplings 


Lepton-flavour violation arises as in Ref. [32], with vks 
playing the role of the non-SM doublet. The Yukawa 
couplings are given by 


^ =-Lf{Cff^3 + Yf-^2)iyR,. 

~ + Yfi'i'2) Ci -b h.C. 


(81) 


and Y® are diagonal due to the 17(1)' symmetry, 
while are given by 

/O 0 0 \ /O 0 0\ 

T'' = 0 0 63 , r = 0 0 0 . (82) 

yo 0 0 / yo Oj 

The right-handed neutrino mass matrix takes again the 
form of Eq. ( |36| ) with two vanishing entries. However, due 
to the non-diagonal Dirac matrices, the active neutrino 
mass matrix no longer features two vanishing minors but 
only one, softening the fine-tuning to obtain valid mix¬ 
ing parameters. We again expect a quasi-degenerate neu¬ 
trino mass spectrum and a close-to-maximal atmospheric 
mixing angle |32j but lose the very specihc predictions we 
obtained in the 2HDM from Sec. III Dl 

Diagonalization of the charged-lepton mass matrix re¬ 
quires only a small 23 rotation of (/r_R,T/j) by an angle 
Or ~ cos f323^Tfiv/^/2mr, while the left-handed angle is 
suppressed by m^/iriT. The LEV coupling of h is then 
approximately given by 


C D 


-Or 


rUr COs(q 23 " 1^23,) 
V COS /323 sin hs 


tPr^H -b h.C. 


(83) 


A non-zero Or hence induces the decay h —>■ pr, while 
any charged LEV involving electrons is forbidden (in the 
limit of zero neutrino masses). The scalars - H 2 , Ag, 
and H 2 - also couple off diagonally to leptons, but their 
effect in LEV observables is small [53] ■ The Z' couplings 
are given by 

C V + (84) 

with ~ diag(0,1, —I) and 

/O 0 0 \ 

Tf 0 1 20 r ]. (85) 

yo 20 r - 1 ) 

Again, a non-zero Or induces LEV in the p-r sector, 
mediated by the Z'. 


C. Lepton-flavour violation 


Concerning LEV we can directly rely on the analysis 
of Refs. [321153]. The branching ratio for h ^ pr reads 


Br [h —)■ pr] ~ 


~ 1 % 



mh Iph |2 

SttTsm 

y / cos(a23 - /323) 

J V 0.2 


/ tan/323 
^ 20 


2 

( 86 ) 


where Tsm — 4.1 MeV is the decay width in the SM for 
a 125 GeV Brout-Englert-Higgs boson and is given 
by the rPRph prefactor in Eq. (831. 

LEV mediated by Z' most importantly induces the de¬ 
cay r —)■ 3/x, with t ^ pj suppressed by an additional 
factor 2aEM/7r [33]. The branching ratio is given by 


Br [t —)■ 3/t] 


mlOl 9 '^ 

1287r3r^ m|, 


10 "® 


/ G.GTeV V 
\ 0 .l) \mz'/g'J 
(87) 













12 


COS(tt23“/^23) ~ 0.25, 3—1/3 
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FIG. 6: Allowed regions in the /g'-sm{6ii) plane for a = 
1/3: the horizontal stripes correspond to h —5 /rr (Icr) for 
tan/323 = 70 , 40 and cos(a 23 — P 23 ) = 0.25, and (light) blue 
stands for (future) r —>■ 3^ limits at 90% C.L. The gray regions 
are excluded by the 2 a range for (see Eq. (561). In this 
range, ATLAS limits constrain mz' > 2.5TeV (see Fig. |^. 


which has to be compared to the current upper limit of 
1.2 X 10“® at 90% C.L. which is obtained from combining 
data from Belle and BaBar [94]. This limit can most 
likely be improved by an order of magnitude to 10“® in 
the future 

In the previous sections, we have seen that a resolution 
of the i?-meson anomalies - indicated through a non-zero 
Cg (Eq. - requires mz' /g' to be in the TeV range 

(Fig. §. In Fig. 1^ we show the exclusion limits from 
T ^ together with the preferred region for h ^ ^lt 
and the Cg constraints on mz'jg' ■ The important part 
is the upper limit on raz' jg' from Cg. With a non-zero 
value for 0^ required by /i —?> //r, we can then predict a 
rate for r —>• 3// mediated by the Z'. For this we express 
mz' Ig' in terms of Cg and Qr in Br[h —>■ p,r] to arrive at 


Br \t 3//] ~ 4.6 x 10"^ 


C| cos^ /323 sin^ /ggs 
a2 cos2(a23 - fe) 


Br[h —)• fj,T]. 


( 88 ) 


We remind the reader that the angles 023 and P 23 do 
not correspond to the 2HDM angles from Sec. |n|but to 
those from Refs. |3H[33]. Using the 2a lower limits on 
Cg (Eq. ( 56 l) and h ^ p,T (Eq. ( pj)), as well as the LHC 
constraint7cos(a23 — /323)| < 0.4 laiTS], we can predict 


Br [r ^ 3/i] > 9.3 x 10"® ( ^ , (89) 

^ ' Vtan/323/ 

working in the large tan/323 limit and setting a = 1/3. 
The current bound is then tan /323 ^ 9, while the future 


reach goes above tan/323 ~ 30. Using the Itr limits for Cg 
and h ^ pLT gives a current (future) bound of 30 (104) 
on tan/323. This is much stronger than the prediction 
of Ref. [33] in a model with vector-like quarks, where 
Icr limits only implied a future reach up to tan/3 ~ 60 
(using the updated value for h ^ gr from Eq. (§). The 
3HDM with gauged horizontal C(l)' charges studied here 
is hence more tightly constrained than the 2HDM with 
vector-like quarks [55] . 

Equation (89) is the main prediction of the simultane¬ 
ous explanation of the H-meson anomalies in connection 
with h —>■ pLT. Note that in addition to the mz'!g' limits 
from Cg, ATLAS constrains mz' vs. g' (Fig. |^. For the 
parameters in Fig. [^ this imposes the additional bound 
mz' > 2.5 TeV (or g' > 0.65), which puts the C(l)' Lan¬ 
dau pole below roughly 3 x 10^^ GeV for a = 1/3. 


V. CONCLUSIONS AND OUTLOOK 

In this paper we proposed a model with multiple 
scalar doublets and a horizontal U{\)' gauge symmetry 
in which all three LHC anomalies in the flavour sector 
{B K*p,'^pL~, R{K) and h —^ //r) can be explained 
simultaneously. Compared to previous explanations, our 
model does not require vector-like quarks charged un¬ 
der the new gauge group. The spontaneously broken 
anomaly-free C/(l)' gauge symmetry is generated by 

Q' = {L^ - Lr) - a(Bi + B 2 - 2 B 3 ), a G Q , (90) 

which leads to successful fermion-mixing patterns. In 
particular, it generates a large (small) atmospheric (re¬ 
actor) mixing angle in the lepton sector and explains the 
almost decoupled third quark generation. The univer¬ 
sal charges the quarks of the first two generations allow 
for the generation of the Cabibbo angle without danger¬ 
ously large effects in Kaon mixing, and the neutralness of 
electrons under the C(l)' symmetry softens constraints 
without fine-tuning. 

The observed quark mixing of the CKM matrix re¬ 
quires the f7(l)' to be broken with a second scalar doublet 
with U{iy charge —a, which leads to flavour-violating 
couplings of the Z' and of the scalars, giving simulta¬ 
neously a natural explanation for the smallness of Vub 
and Vcb- Scalar contributions to Bg-Bg mixing typi¬ 
cally require a — /3 ~ 7r/2, which is, however, relaxed 
for rriA < rnn- The anomalies in H —>■ and 

R{K) can be explained with a TeV-scale Z' boson and 
a < 1 while satisfying Hg-Hs-mixing constraints and lim¬ 
its from direct Z' searches at the LHC. Future LHC and 
FCC (Future Circular Collider) searches are very inter¬ 
esting for our model as they might strengthen the current 
limits or lead to the discovery of the Z' boson. 

Introducing a third scalar doublet, with 11(1)' charge 
—2, gives rise to the decay h ^ ptr in complete analogy to 
Refs. 1321133]. Together with the large Z' effect necessary 
to resolve B —)• K*//+//“ and R{K), the decay h ^ pLT 
then allows us to predict a rate for r —)■ 3//, depending 
on tan /3 and cos(a — /3), potentially measurable in future 
experiments. 
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Note added: During the publication process of this ar¬ 
ticle, new LHCb results were presented at the Rencon¬ 
tres de Moriond Electroweak Session 2015 which hint at 
a confirmation of the anomaly in H —>■ K* . The 
global fit now prefers new physics in Cg ^ over the Stan¬ 
dard Model by 4.3 (t [HH] , 


Appendix A: Other horizontal symmetries 

Demanding a universal U{1)' quark coupling to the 
hrst two generations and a good flavor symmetry in the 
lepton sector does not uniquely single out our model with 
Bi -\- B 2 — 2 B 3 and — Lt-. It is well known that, be¬ 


sides — Lt (which is connected to quasi-degenerate 
neutrinos), the symmetries Lg and Lg — — Lt are good 

zeroth-order approximations for a neutrino mass matrix 
with normal and inverted hierarchy, respectively |55j . 
Since these two are anomalous, one can consider B — 
3Le EZl EH] or B -\- 3{Lg — — Lt) HOO) as well- 

motivated gauge symmetries. With non-universal quark 
charges - but universal in the first two generations - we 
can consider B 3 — Lg or B 3 (Lg — — Lt) as anomaly- 

free gauge symmetries that provide a successful neutrino 
mixing, single out the third quark generation, and lead to 
LFV in both lepton and quark sectors. Quite analogously 
one can consider the lepton symmetries from Ref. |59j 
with a non-universal quark charge, which lead to predic¬ 
tive texture zeros and vanishing minors in the neutrino 
mass matrix; for example, 3 R 3 -\- Lg — 3L^ — Lt gener¬ 
ates the texture zeros (TOjy)ii = 0 = after seesaw. 

While all these symmetries are interesting in their own 
right, they are not useful for our purpose because the 
Z' coupling to quarks compared to muons is rather large 
(and not adjustable), so it becomes difficult to generate 
a large Cg without violating L-L-mixing bounds. In ad¬ 
dition, any Z' that couples to electrons unavoidably suf¬ 
fers from stringent LEP constraints m- The horizontal 
gauge symmetry (L^ — Lt) — a{Bi +B 2 — 2 B 3 ) chosen in 
this paper is hence a remarkably good choice to address 
the existing hints for flavour violation. 
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